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System, Method, and Computer Program Product for the Visualization 
and Interactive Processing and Analysis of Chemical Data 

Background of the Invention 

Field of the Invention 

The present invention is generally directed to displaying and processing data using a 
computer, and more particularly directed to visualizing and interactively processing chemical 
compounds using a computer. 

Related Art 

Currently, research to identify chemical compounds with useful properties (such as 
paints, finishes, plasticizers, surfactants, scents, drugs, herbicides, pesticides, veterinary 
products, etc.) often includes the synthesis/acquisition and analysis of large libraries of 
chemical compounds. More and more, combinatorial chemical libraries are being 
synthesized/acquired and analyzed to conduct this research. 

A combinatorial chemical library is a collection of diverse chemical compounds 
generated by either chemical synthesis or biological synthesis by combining a number of 
chemical "building blocks" such as reagents. For example, a linear combinatorial chemical 
library such as a polypeptide library is formed by combining a set of chemical building 
blocks called amino acids in every possible way for a given compound length (i.e., the 
number of amino acids in a polypeptide compound). Millions of chemical compounds 
theoretically can be synthesized through such combinatorial mixing of chemical building 
blocks. For example, one commentator has observed that the systematic, combinatorial 
mixing of 100 interchangeable chemical building blocks results in the theoretical synthesis of 
100 million tetrameric compounds or 10 billion pentameric compounds (Gallop et a/., 
"Applications of Combinatorial Technologies to Drug Discovery, Background and Peptide 
Combinatorial Libraries," J. Med. Chem. 37, 1233-1250 (1994)). 
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Advanced research in this area often involves the use of directed diversity libraries. 
A directed diversity library is a large collection of chemical compounds having 
properties/features/characteristics that match some prescribed properties. The generation, 
analysis, and processing of directed diversity libraries are described in U.S. Patent Nos. 
5 5,463,564; 5,574,656; and 5,684,711, and pending U.S. Application titled "SYSTEM, 
METHOD AND COMPUTER PROGRAM PRODUCT FOR IDENTIFYING CHEMICAL 
COMPOUNDS HAVING DESIRED PROPERTIES," Atty. Docket No. 1503.0200001, all of 
which are herein incorporated by reference in their entireties. 

In conducting such research, it would be very valuable to be able to compare the 
10 properties, features, and other identifying characteristics of compounds. For example, 
suppose that a researcher has identified a compound X that exhibits some useful properties. 
It would aid the researcher greatly if he could identify similar compounds, since those similar 
compounds might also exhibit those same useful properties. 

It would also help a researcher in his work to be able to easily synthesize 
15 compounds, or retrieve compounds from a chemical inventory. Further, it would greatly aid 
a researcher to be able to interactively analyze and otherwise process chemical compounds. 

Summary of the Invention 

Briefly stated, the present invention is directed to a system, method, and computer 
20 program product for visualizing and interactively analyzing data relating to chemical 
compounds. The invention operates as follows. A user selects a plurality of compounds to 
map, and also selects a method for evaluating similarity/dissimilarity between the selected 
compounds. A non-linear map is generated in accordance with the selected compounds and 
the selected method. The non-linear map has a point for each of the selected compounds, 
25 wherein a distance between any two points is representative of similarity/dissimilarity 
between the corresponding compounds. A portion of the non- linear map is then displayed. 
Users are enabled to interactively analyze compounds represented in the non- linear map. 

Further features and advantages of the present invention, as well as the structure and 
operation of various embodiments of the present invention, are described in detail below with 
30 reference to the accompanying drawings. In the drawings, like reference numbers indicate 
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identical or functionally similar elements. Also, the leftmost digit(s) of the reference 
numbers identify the drawings in which the associated elements are first introduced. 

Brief Description of the Figures 

The file of this patent contains at least one drawing executed in color. Copies of this 
patent with color drawing(s) will be provided by the Patent and Trademark Office upon 
request and payment of the necessary fee. 

The present invention will be described with reference to the accompanying 
drawings, wherein: 

FIG. 1 illustrates a block diagram of a computing environment according to an 
embodiment of the invention; 

FIG. 2 is a block diagram of a computer useful for implementing components of the 
invention; 

FIG. 3 is a flowchart representing the operation of the invention in visualizing and 
interactively processing non-linear maps according to an embodiment of the invention; 

FIG. 4 is a flowchart representing the manner in which a non-linear map is 
generated according to an embodiment of the invention; 

FIG. 5 illustrates a structure browser window according to an embodiment of the 
invention; 

FIG. 6 illustrates a compound visualization non-linear map window according to an 
embodiment of the invention; 

FIG. 7 is used to describe a zoom function of the present invention; 

FIG. 8 illustrates a dialog used to adjust properties of a set containing one or more 
compounds; 

FIGS. 9 and 10 are used to describe the compound visualization non-linear map 
window according to an embodiment of the invention; 

FIG. 1 1 is a flowchart illustrating the operation of the invention where a compound 
visualization non-linear map window is used as a source in an interactive operation; 

FIG. 12 is a flowchart illustrating the operation of the invention where a compound 
visualization non-linear map window is used as a target in an interactive operation; 
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F1G. 13 conceptually illustrates an interactive operation where a compound 
visualization non-linear map window is used as a source; and 

FIG. 14 conceptually illustrates an interactive operation where a compound 
visualization non-linear map window is used as a target. 
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Detailed Description of the Preferred Embodiments 

Table of Contents 

1 . Overview of the Present Invention 

2. Structure of the Invention 

3. Implementation Embodiment of the Invention 

4. Overview of Multidimensional Scaling (MDS) and Non-Linear Mapping (NLM) 

4. 1 Procedure Suitable for Relatively Small Data Sets 

4.2 Procedure Suitable for Large Data Sets 

5. Evaluation Properties (Features) and Distance Measures 

5. 1 Evaluation Properties Having Continuous or Discrete Real Values 

5.2 Distance Measure Where Values of Evaluation Properties Are Continuous or 
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5.3 Evaluation Properties Having Binary Values 

5.4 Distance Measures Where Values of Evaluation Properties Are Binary 

6. Scaling of Evaluation Properties 
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10. Examples 
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7. Overview of the Present Invention 

The present invention is directed to a computer-based system, method, and/or 
computer program product for visualizing and analyzing chemical data using interactive 
multi-dimensional (such as 2- and/or 3-dimensional) non-linear maps. In particular, the 
5 invention employs a suite of non-linear mapping algorithms to represent chemical 
compounds as objects in preferably 2D or 3D Euclidean space. 

According to the invention, the distances between objects in that space represent the 
similarities and/or dissimilarities of the corresponding compounds (relative to selected 
properties or features of the compounds) computed by some prescribed method. The 
10 resulting maps are displayed on a suitable graphics device (such as a graphics terminal, for 
example), and interactively analyzed to reveal relationships between the data, and to initiate 
an array of tasks related to these compounds. 

2. Structure of the Invention 

FIG. 1 is a block diagram of a computing environment 102 according to a preferred 
1 5 embodiment of the present invention. 

A chemical data visualization and interactive analysis module 104 includes a map 
generating module 106 and user interface modules 108. The map generating module 106 
determines distances between chemical compounds relative to one or more selected 
properties or features (herein sometimes called evaluation properties or features) of the 
20 compounds. The map generating module 106 performs this function by retrieving and 
analyzing data on chemical compounds and reagents from reagent and compound databases 
122. These reagent and compound databases 122 store information on chemical compounds 
and reagents of interest. 

The reagent and compound databases 122 are part of databases 120, which 
25 communicate with the chemical data visualization and interactive analysis module 104 via a 
communication medium 1 1 8. The communication medium 1 1 8 is preferably any type of data 
communication means, such as a data bus, a computer network, etc. 

The user interface modules 108, which include a map viewer 1 12 and optionally a 
structure browser 1 10, displays a preferably 2D or 3D non-linear map on a suitable graphics 
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device. The non-linear map includes objects that represent the chemical compounds, where 
the distances between the objects in the non-linear map are those distances determined by the 
map generating module 106. The user interface modules 108 enable human operators to 
interactively analyze and process the information in the non-linear map so as to reveal 
relationships between the data, and to initiate an array of tasks related to the corresponding 
compounds. 

The user interface modules 108 enable users to organize compounds as collections 
(representing, for example, a combinatorial library). Information pertaining to compound 
collections are preferably stored in a collection database 124. Information on reagents that 
are mixed to form compound collections are preferably stored in a library database 126. 

Input Device(s) 114 receive input (such as data, commands, queries, etc.) from 
human operators and forward such input to, for example, the chemical data visualization and 
interactive analysis module 104 via the communication medium 118. Any well known, 
suitable input device can be used in the present invention, such as a keyboard, pointing 
device (mouse, roller ball, track ball, light pen, etc.), touch screen, voice recognition, etc. 
User input can also be stored and then retrieved, as appropriate, from data/command files. 

Output Device(s) 116 output information to human operators. Any well known, 
suitable output device can be used in the present invention, such as a monitor, a printer, a 
floppy disk drive or other storage device, a text-to-speech synthesizer, etc. 

As described below, the present invention enables the chemical data visualization 
and interactive analysis module 104 to interact with a number of other modules, including 
but not limited to one or more map viewers 112, NMR (nuclear magnetic resonance) 
widget/module 130, structure viewers 110, MS (mass spectrometry) widget/module 134, 
spreadsheets 136, QSAR (Quantitative Structure-Activity Relationships) module 138, an 
experiment planner 140, property prediction programs 142, active site docker 144, etc. 
These modules communicate with the chemical data visualization and interactive analysis 
module 104 via the communication medium 118. 
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J. Implementation Embodiment of the Invention 

Components shown in the computing environment 102 of FIG. 1 (such as the 
chemical data visualization and interactive analysis module 104) can be implemented using 
one or more computers, such as an example computer 202 shown in FIG. 2. 
5 The computer 202 includes one or more processors, such as processor 204. 

Processor 204 is connected to a communication bus 206. Various software embodiments are 
described in terms of this example computer system. After reading this description, it will 
become apparent to a person skilled in the relevant art(s) how to implement the invention 
using other computer systems and/or computer architectures. 

10 Computer 202 also includes a main memory 208, preferably random access memory 

(RAM), and can also include one or more secondary storage devices 210. Secondary storage 
devices 210 can include, for example, a hard disk drive 212 and/or a removable storage drive 
214, representing a floppy disk drive, a magnetic tape drive, an optical disk drive, etc. 
Removable storage drive 214 reads from and/or writes to a removable storage unit 216 in a 

15 well known manner. Removable storage unit 216 represents a floppy disk, magnetic tape, 
optical disk, etc. which is read by and written to by removable storage drive 214. Removable 
storage unit 216 includes a computer usable storage medium having stored therein computer 
software and/or data. 

In alternative embodiments, the computer 202 can include other similar means for 
20 allowing computer programs or other instructions to be loaded into computer 202. Such 
means can include, for example, a removable storage unit 220 and an interface 218. 
Examples of such can include a program cartridge and cartridge interface (such as that found 
in video game devices), a removable memory chip (such as an EPROM, or PROM) and 
associated socket, and other removable storage units 220 and interfaces 218 which allow 
25 software and data to be transferred from the removable storage unit 220 to computer 202. 

The computer 202 can also include a communications interface 222. 
Communications interface 222 allows software and data to be transferred between computer 
202 and external devices. Examples of communications interface 222 include, but are not 
limited to a modem, a network interface (such as an Ethernet card), a communications port, a 
30 PCMCIA slot and card, etc. Software and data transferred via communications interface 222 
are in the form of signals which can be electronic, electromagnetic, optical or other signals 
capable of being received by communications interfac 222. 
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In this document, the term "computer program product" is used to generally refer to 
media such as removable storage units 216, 220, a hard drive 212 that can be removed from 
the computer 202, and signals carrying software received by the communications interface 
222. These computer program products are means for providing software to the computer 

5 202. 

Computer programs (also called computer control logic) are stored in main memory 
and/or secondary storage devices 210. Computer programs can also be received via 
communications interface 222. Such computer programs, when executed, enable the 
computer 202 to perform the features of the present invention as discussed herein. In 
1 0 particular, the computer programs, when executed, enable the processor 204 to perform the 
features of the present invention. Accordingly, such computer programs represent controllers 
of the computer 202. 

In an embodiment where the invention is implemented using software, the software 
can be stored in a computer program product and loaded into computer 202 using removable 

15 storage drive 214, hard drive 212, and/or communications interface 222. The control logic 
(software), when executed by the processor 204, causes the processor 204 to perform the 
functions of the invention as described herein. 

In another embodiment, the automated portion of the invention is implemented 
primarily in hardware using, for example, hardware components such as application specific 

20 integrated circuits (ASICs). Implementation of the hardware state machine so as to perform 
the functions described herein will be apparent to persons skilled in the relevant art(s). 

In yet another embodiment, the invention is implemented using a combination of 
both hardware and software. 

The computer 202 can be any suitable computer, such as a computer system running 

25 an operating system supporting a graphical user interface and a windowing environment. A 
suitable computer system is a Silicon Graphics, Inc. (SGI) workstation/server, a Sun 
workstation/server, a DEC workstation/server, an IBM workstation/server, an IBM 
compatible PC, an Apple Macintosh, or any other suitable computer system, such as one 
using one or more processors from the Intel Pentium family, such as Pentium Pro or Pentium 

30 II. Suitable operating systems include, but are not limited to, IRIX, OS/Solaris, Digital Unix, 
AIX, Microsoft Windows 95/NT, Apple Mac OS, or any other operating system supporting a 
graphical user interface and a windowing environment. For example, in a preferred 
embodim nt the program may be 



SUBSTITUTE SHEET (RULE 26) 



WO 98/20459 



PCT/US97/20919 



-10- 

implemented and run on an Silicon Graphics Octane workstation running the IRIX 6.4 
operating system, and using the Motif graphical user interface based on the X Window 
System. 

4. Overview of Multidimensional Scaling (MDS) and Non-Linear Mapping (NLM) 

According to the present invention, multidimensional scaling (MDS) and non-linear 
mapping (NLM) techniques are used to generate the non-linear map (i.e., the non-linear 
map) that includes objects, where the objects represent chemical compounds, and the 
distances between the objects are indicative of the similarities and dissimilarities between the 
corresponding compounds. MDS and NLM are described in this section. 

MDS and NLM were introduced by Torgerson, Phychometrika, 17:401 (1952); 
Kniskai, Psychometrika, 29:115 (1964); and Sammon, IEEE Trans, Compute C- 18:40 1 
(1969) as a means to generate low-dimensional representations of psychological data. 
Multidimensional scaling and non-linear mapping are reviewed in Schiffman, Reynolds and 
Young, Introduction to Multidimensional Scaling, Academic Press, New York (1981); 
Young and Hamer, Multidimensional Scaling: History, Theory and Applications, Erlbaum 
Associates, Inc., Hillsdale, NJ (1987); and Cox and Cox, Multidimensional Scaling, Number 
59 in Monographs in Statistics and Applied Probability, Chapman-Hall (1994). The contents 
of these publications are incorporated herein by reference in their entireties. 

4 J Procedure Suitable for Relatively Small Data Sets 

MDS and NLM (these are generally the same, and are hereafter collectively referred 
to as MDS) represent a collection of methods for visualizing proximity relations of objects 
by distances of points in a low-dimensional Euclidean space. Proximity measures are 
reviewed in Hartigan, J. Am. Statist Ass., 62:1 140 (1967), which is incorporated herein by 
reference in its entirety. In particular, given a finite set of vectorial or other samples A = i 
= 1, k}, a distance function dy = d^, with a* a^ € A, which measures the similarity 
and dissimilarity between the i-th and j-th objects in A, and a set of images X = {x if X*; x t 
e of A on an m-dimensional display plane (3* m being an m dimensional vector of real 
numbers), the objective is to place * onto the display plane in such a way that their Euclidean 
distances ||x, - Xj|| 
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approximate as closely as possible the corresponding values d^. This projection, which can 
only be made approximately, is carried out in an iterative fashion by minimizing an error 
function which measures the difference between the distance matrices of the original and 
projected vector sets. Several such error functions have been proposed, most of which are of 
the least-squares type, including Kjuskal's 'stress*: 



' v it ir 



\ 



EQ. 1 



Sammon's error criterion: 




E', 

EQ. 2 



and Lingoes' alienation coefficient: 



10 



K = 



N 



E« 9 



EQ.3 



where 6^ = ||Xj - Xj|| is the Euclidean distance between the images and x } on the display 
plane. Generally, the solution is found in an iterative fashion by (I) computing or retrieving 
from a database the distances d {i ; (2) initializing the images iji (3) computing the distances of 
the images 
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6 and the value of the error function (e.g. S, E or K in EQ. 1-3 above); (4) computing a new 
configuration of the images x 5 using a gradient descent procedure, such as Kruskal's linear 
regression or Guttmans rank-image permutation; and (5) repeating steps 3 and 4 until the 
error is minimized within some prescribed tolerance. 
5 For example, the Sammon algorithm minimizes EQ. 2 by iteratively updating the 

coordinates x, using Eq 4: 

x (m + l)=x (m)-AA (m) A 

pq\ * pq\ I pq\ ' £Q 4 



where m is the iteration number, is the q-th coordinate of the p-th image Xp, X is the 
learning rate, and 
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The partial derivatives in EQ. 5 are given by: 
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The non-linear mapping is obtained by repeated evaluation of EQ. 2, followed by 
modification of the coordinates using EQ. 4 and 5, until the error is minimized within a 
prescribed tolerance. 

4.2 Procedure Suitable for Large Data Sets 

The general refinement paradigm described in Section 4.1 is suitable for relatively 
small data sets, but has one important limitation that renders it impractical for large data sets. 
This limitation stems from the fact that the computational effort required to compute the 
gradients scales to the square of the size of the data set For relatively large data sets, this 
quadratic time complexity makes even a partial refinement intractable. 

According to the present invention, the following approach is used for large data sets. 
This approach is to use iterative refinement based on 'instantaneous' errors. As in the 
approach described in Section 4.1, this approach of Section 4.2 starts with an initial 
configuration of points generated at random or by some other procedure (as described below 
in Section 7). This initial configuration is then continuously refined by repeatedly selecting 
two points i, j, at random, and modifying their coordinates on the non-linear map according 
to Eq. 8: 

x t (t + 1 ) =A^(0^ ; (0,4,) EQ 8 



where t is the current iteration, Xj(t) and Xj(t) are the current coordinates of the i-th and j-th 
points on the non-linear map, Xj(t+1) are the new coordinates of the i-th point on the non- 
linear map, and is the true distance between the i-th and j-th points that we attempt to 
approximate on the non-linear map (see above). J{.) in EQ. 8 above can assume any 
functional form. Ideally, this function should try to minimize the difference between the 
actual and target distance between the i-th and j-th points. For example,/.) may be given by 
EQ. 9: 



(d -6 (0), 



EQ.9 
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where t is the iteration number, 6 V) = ||Xj(t) - x } (t)\\ , and X(x) is an adjustable parameter, 
referred to hereafter as the 'learning rate.' 

An analogous equation has been suggested by Kohonen for the training of self- 
organizing maps (Kohonen, Self-Organizing Maps, Springer- Verlag, Berlin (1995)), 
5 incorporated herein by reference in its entirety. This process is repeated for a fixed number 
of cycles, or until some global error criterion is minimized within some prescribed tolerance. 
A large number of iterations are typically required to achieve statistical accuracy. 

The method described above is generally reminiscent of Kohonen's self-organizing 
principle (Kohonen, Biological Cybernetics, 43:59 (1982)) and neural network back- 

10 propagation training (Werbos, Beyond Regression: New Tools for Prediction and Analysis in 
the Behavioral Sciences, PhD Thesis^ Harvard University, Cambridge, MA (1974)), and 
Rumeihart and McClelland, Eds., Parallel Distributed Processing: Explorations in the 
Microstructure of Cognition, Vol. 1, MIT Press, Cambridge, MA (1986)), all of which are 
incorporated herein by reference in their entireties. 

15 The learning rate X(\) in EQ. 9 plays a key role in ensuring convergence. If X is too 

small, the coordinate updates are small, and convergence is slow. If, on the other hand, X is 
too large, the rate of learning may be accelerated, but the non-linear map may become 
unstable (i.e. oscillatory). Typically, X ranges in the interval [0, 1] and may be fixed, or it 
may decrease monotonically during the refinement process. Moreover, X may also be a 

20 function of i, j and/or d ijf and can be used to apply different weights to certain objects, 
distances and/or distance pairs. For example, X may be computed by EQ. 10: 



X . -X i 
EQ. 10 



orEQ. 11: 
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EQ. 1 1 

where X mM and k min are the (unweighted) starting and ending learning rates such that k mw * min 
e [0,1], T is the total number of refinement steps (iterations), t is the current iteration 
number, and a is a constant scaling factor. EQ. 10 and It have the effect of decreasing the 
correction at large separations, thus creating a non-linear map which preserves short-range 
interactions more faithfully than long-range ones. Weighting is discussed in greater detail 
below. Because of the general resemblance of the training process described above to 
Kohonen's self-organizing principle, these maps shall sometimes be herein called 'Self- 
Organizing Non-Linear Maps.' 

One of the main advantages of this approach is that it makes partial refinements 
possible. It is often sufficient that the pair-wise dissimilarities are represented only 
approximately to reveal the general structure and topology of the data. Unlike traditional 
MDS, this approach allows very fine control of the refinement process. Moreover, as the 
non-linear map self-organizes, the pair-wise refinements become cooperative, which partially 
alleviates the quadratic nature of the problem. 

The general usefulness of multi-dimensional scaling stems from the fact that data in 9t* are 
almost never d-dimensional. Although scaling becomes more problematic as the true 
dimensionality of the space increases, the presence of structure in the data is very frequently 
reflected on the resulting map. Of course, one can easily conceive of situations where MDS 
is not effective, particularly when the data is random and truly hyper-dimensional. 
Fortunately, these situations rarely arise in practice, as some form of structure is always 
present in the data, particularly data related to molecular structure and function. 

The embedding procedure described above does not guarantee convergence to the 
global minimum (i.e., the most faithful embedding in a least-squares sense). If so desired, the 
refinement process may be repeated a number of times from different starting configurations 
and/or random number seeds. It should also be pointed out that the absolute coordinates in 
the non-linear map carry no physical significance. What is important are the relative 
distances 
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between points, and the general structure and topology of the data (presence, density and 
separation of clusters, etc.). 

The method described above is ideally suited for both metric and non-metric scaling. 
The latter is particularly useful when the (dis)similarity measure is not a true metric, i.e. it 
5 does not obey the distance postulates and, in particular, the triangle inequality (such as the 
Tanimoto coefficient, for example). Although an "exact* projection is only possible when the 
distance matrix is positive definite, meaningful projections can still be obtained even when 
this criterion is not satisfied. As mentioned above, the overall quality of the projection is 
determined by a sum-of-squares error function such as those shown in EQ. 1-3. 

1 0 5. Evaluation Properties (Features) and Distance Measures 

As mentioned above, the distances d Sj between chemical compounds are computed 
according to some prescribed measure of molecular 'similarity'. This similarity can be based 
on any combination of properties or features of the compounds. For example, the similarity 
measure may be based on structural similarity, chemical similarity, physical similarity, 
15 biological similarity, and/or some other type of similarity measure which can be derived 
from the structure or identity of the compounds. Under the system of the present invention, 
any similarity measure can be used to construct the non-linear map. The properties or 
features that are being used to evaluate similarity or dissimilarity among compounds are 
sometimes herein collectively called "evaluation properties." 

20 5.7 Evaluation Properties Having Continuous or Discrete Real Values 

As noted above, in a preferred embodiment of the present invention, the similarity 
measure may be derived from a list of physical, chemical and/or biological properties (i.e., 
evaluation properties) associated with a set of compounds. Under this formalism, the 
compounds are represented as vectors in multi-variate property space, and their similarity 
25 may be computed by some geometrical distance measure. 

In a preferred embodiment, the property space is defined using one or more 
molecular features (descriptors). Such molecular features may include topological indices, 
physicochemical 
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properties, electrostatic field parameters, volume and surface parameters, etc. For example, 
these features may include, but are not limited to, molecular volume and surface areas, dipole 
moments, octanol-water partition coefficients, molar refractivities, heats of formation, total 
energies, ionization potentials, molecular connectivity indices, 2D and 3D auto-correlation 

5 vectors, 3D structural and/or pharmacophoric parameters, electronic fields, etc. However, it 
should be understood that the present invention is not limited to this embodiment. For 
example, molecular features may include the observed biological activities of a set of 
compounds against an array of biological targets such as enzymes or receptors (also known 
as affinity fingerprints). In fact, any vectorial representation of chemical data can be used in 

1 0 the present invention. 

5.2 Distance Measure Where Values of Evaluation Properties Are Continuous 
or Discrete Real Numbers 



A "distance measure" is some algorithm or technique used to determine the 
difference between compounds based on the selected evaluation properties. The particular 
15 distance measure that is used in any given situation depends, at least in part, on the set of 
values that the evaluation properties can take. 

For example, where the evaluation properties can take real numbers as values, then a 
suitable distance measure is the Minkowski metric, shown in EQ. 12: 



EQ. 12 



20 where k is used to index the elements of the property vector, and re[ 1 , »). For r = 1 .0, EQ. 1 2 
is the city-block or Manhattan metric. For r = 2.0, EQ. 12 is the ordinary Euclidean metric. 
For r = «, EQ. 12 is the maximum of the absolute coordinate distances, also referred to as the 
'dominance' metric, the 'sup' metric, or the 'ultrametric' distance. For any value of re[I, ~), 
it can be shown that the Minkowski metric is a true metric, i.e. it obeys the distance 

25 postulates and, in particular, the triangle inequality. 
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5.3 Evaluation Properties Having Binary Values 

Alternatively, the evaluation properties of the compounds may be represented in a 
binary form (i.e., either a compound has or does not have an evaluation property), where 
each bit is used to indicate the presence or absence (or potential presence or absence) of 
5 some molecular feature or characteristic. For example, compounds may be encoded using 
substructure keys where each bit is used to denote the presence or absence of a specific 
structural feature or pattern in the target molecule. Such features include, but are not limited 
to, the presence, absence or minimum number of occurrences of a particular element (e.g. the 
presence of at least I, 2 or 3 nitrogen atoms), unusual or important electronic configurations 

10 and atom types (e.g. doubly-bonded nitrogen or aromatic carbon), common functional groups 
such as alcohols, amines etc, certain primitive and composite rings, a pair or triplet of 
pharmacophoric groups at a particular separation in 3-dimensional space, and 'disjunctions' 
of unusual features that are rare enough not to worth an individual bit, yet extremely 
important when they do occur (typically, these unusual features are assigned a common bit 

1 5 that is set if any one of the patterns is present in the target molecule). 

Alternatively, the evaluation properties of compounds may be encoded in the form of 
binary fingerprints, which do not depend on a predefined fragment or feature dictionary to 
perform the bit assignment. Instead, every pattern in the molecule up to a predefined limit is 
systematically enumerated, and serves as input to a hashing algorithm that turns W a small 

20 number of bits at pseudo-random positions in the bitmap. Although it is conceivable that two 
different molecules may have exactly the same fingerprint, the probability of this happening 
is extremely small for all but the simplest cases. Experience suggests that these fingerprints 
contain sufficient information about the molecular structures to permit meaningful similarity 
comparisons. 

25 5.4 Distance Measures Where Values of Evaluation Properties Are Binary 

A number of similarity (distance) measures can be used with binary descriptors (i.e., 
where evaluation properties are binary or binary fingerprints). The most frequently used ones 
are the normalized Hamming distance: 
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!I _ \XOR(x t y)\ 

. N E0..3. 

which measures the number of bits that are different between x and y, the Tanimoto or 
Jaccard coefficient: 



T _ [iND(x,y)\ 

^ OR ^\ EQ. .4 

which is a measure of the number of substructures shared by two molecules relative to the 
ones they could have in common, and the Dice coefficient: 

2[iND(x,y)\ 

w ^ EQ. 15 

In the equations listed above, AND(x, y) is the intersection of binary sets x and y (bits that 
are 'on' in both sets), IOR(x, y) is the union or 'inclusive or* of x and y (bits that are 'on' in 
either x or y), XOR is the 'exclusive or* of x and y (bits that are 'on' in either x or y, but not 
both), |x| is the number of bits that are 'on' in x, and N is the length of the binary sets 
measured in bits (a constant). 

Another popular metric is the Euclidean distance which, in the case of binary sets, 
can be recast in the form: 



E=y/N-\XOR(x t NOT(y))\ fi 



where NOT(y) denotes the binary complement of y. The expression |XOR(x, NOT(y))| 
represents the number of bits that are identical in x and y (either 1 's or O's). The Euclidean 
distance is a good measure of similarity when the binary sets are relatively rich, and is 
mostly used in situations in which similarity is measured in a relative sense. 
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In the examples described above, the distance between two compounds is determined 
using a binary or multivariate representation. However, the system of the present invention is 
not limited to this embodiment. For example, the similarity between two compounds may be 
determined by comparing the shapes of the molecules using a suitable 3-dimensional 
5 alignment method, or it may be inferred by a similarity model defined according to a 
prescribed procedure. For example, one such similarity model may be a neural network 
trained to predict a similarity coefficient given a suitably encoded pair of compounds. Such a 
neural network may be trained using a training set of structure pairs and a known similarity 
coefficient for each such pair, as determined by user input, for example. 

1 0 6. Scaling of Evaluation Properties 

Referring back to EQ. 12, according to the present invention, the features (i.e., 
evaluation properties) may be scaled differently to reflect their relative importance in 
assessing the proximity between two compounds. For example, suppose the user has selected 
two evaluation properties, Property A and Property B. If Property A has a weight of 2, and 
15 Property B has a weight of 10, then Property B will have five times the impact on the 
distance calculation than Property A. 

According to this embodiment of the invention, EQ. 12 may be replaced by EQ. 17: 

V * ' EQ. 17 

where w k is the weight of the k-th property. An example of such a weighting factor is a 
20 normalization coefficient. However, other weighting schemes may also be used. 

According to the present invention, the scaling (weights) need not be uniform 
throughout the entire map, i.e. the resulting map need not be isomorphic. Hereafter, maps 
derived from uniform weights shall be referred to as globally weighted (isomorphic), 
whereas maps derived from non-uniform weights shall be referred to as locally weighted 
25 (non-isomorphic). On locally-weighted maps, the distances on the non-linear map reflect a 
local measure of similarity. That is, 
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what determines similarity in one domain of the non-linear map is not necessarily the same 
with what determines similarity on another domain of the non-linear map. For example, 
locally-weighted maps may be used to reflect similarities derived from a locally- weighted 
case-based learning algorithm. Locally- weighted learning uses locally weighted training to 
5 average, interpolate between, extrapolate from, or otherwise combine training data. Most 
learning methods (also referred to as modeling or prediction methods) construct a single 
model to fit all the training data. Local models, on the other hand, attempt to fit the training 
data in a local region around the location of the query. Examples of local models include 
nearest neighbors, weighted average, and locally weighted regression. Locally- weighted 

10 learning is reviewed in Vapnik, in Advances in Neural Information Processing Systems , 
4:831, Morgan-Kaufman, San Mateo, CA (1982); Bottou and Vapnik, Neural Computation, 
4(6):888 (1992); and Vapnik and Bottou, Neural Computation, 5(6):893 (1993), all of which 
are incorporated herein by reference in their entireties. 

According to the present invention, it is also possible to construct a non-linear map 

1 5 from a distance matrix which is not strictly symmetric, i.e. a distance matrix where d^ * d^. A 
potential use of this approach is in situations where the distance function is defined locally, 
e.g. in a locally weighted model using a point-based local distance function. In this 
embodiment, each training case has associated with it a distance function and the values of 
the corresponding parameters. Preferably, to construct a non- linear map which reflects these 

20 local distance relationships, the distance between two points is evaluated twice, using the 
local distance functions of the respective points. The resulting distances are averaged, and are 
used as input in the non-linear mapping algorithm described above. If the point-based local 
distance functions vary in some continuous or semi-continuous fashion throughout the 
feature space, this approach could potentially lead to a meaningful projection. 

25 7. Improvements to Map Generation Process 

This section describes improvements to the chemical visualization map generation 
process described above. Each of the enhancements described below is under the control of 
the user. That is, the user can elect to perform or not perform each of the enhancements 
discussed below. Alternatively, the invention can be defined so that the below enhancements 
30 are 
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automatically performed, unless specifically overrided by the user (or in some embodiments, 
the user may not have the option of overriding one or more of the below enhancements). 

7.7 Pre~Ordering 

In many cases, the approach described above for generating the non-linear map may 
5 be accelerated by pre-ordering the data using a suitable statistical method. For example, if the 
data is available in vectorial or binary form, the initial configuration of the points on the non- 
linear map may be computed using Principal Component Analysis. In a preferred 
embodiment, the initial configuration may be constructed from the first 3 principal 
components of the feature matrix (i.e. the 3 latent variables which account for most of the 

10 variance in the data). In practice, this technique can have profound effects in the speed of 
refinement. Indeed, if a random initial configuration is used, a significant portion of the 
training time is spent establishing the general structure and topology of the non-linear map, 
which is typically characterized by large rearrangements. If, on the other hand, the input 
configuration is partially ordered, the error criterion can be reduced relatively rapidly to an 

1 5 acceptable level. 

7.2 Localized Refinement 

If the data is highly clustered, by virtue of the sampling process low-density areas 
may be refined less effectively than high-density areas. In a preferred embodiment, this 
tendency may be partially compensated by a modification to the original algorithm which 

20 increases the sampling probability in low-density areas. In one embodiment, the center of 
mass of the non-linear map is identified, and concentric shells centered at that point are 
constructed. A series of regular refinement iterations are then carried out, each time selecting 
points from within or between these shells. This process is repeated for a prescribed number 
of cycles. This phase is then followed by a phase of regular refinement using global 

25 sampling, and the process is repeated. 

As mentioned above, the basic algorithm does not distinguish short- from long-range 
distances. EQ. 10 and 1 1 describe a method to ensure that short-range distances are preserved 
more faithfully than long-range ones through the use of weighting. An alternative (and 
complementary) approach is to ensure that points at close separation are sampled more 
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extensively than points at long separation. A preferred embodiment is to use an alternating 
sequence of global and local refinement cycles, similar to the one described above. In this 
embodiment, a phase of global refinement is initially carried out. At the end of this phase, the 
resulting non-linear map is partitioned into a regular grid, and the points (objects) in each cell 
5 are subjected to a phase of local refinement (i.e. only points from within the same cell are 
compared and refined). Preferably, the number of sampling steps in each cell should be 
proportional to the number of points contained in that cell. This process is highly 
parallelizable. This local refinement phase is then followed by another global refinement 
phase, and the process is repeated for a prescribed number of cycles, or until the embedding 
10 error is minimized within a prescribed tolerance. Alternatively, the grid method may be 
replaced by another suitable method for identifying proximal points, such as a k-d tree, for 
example. 

73 Incremental Refinement 

The approach and techniques described herein may be used for incremental 
1 5 refinement of a map. That is, starting from an organized non-linear map of a set of objects or 
points (compounds), a new set of objects (compounds) may be added without modification of 
the original map. Strictly speaking, this is statistically acceptable if the new set of objects is 
significantly smaller than the original set. In a preferred embodiment, the new set of objects 
may be 'diffused' into the existing map, using a modification of the algorithm described 
20 above. In particular, EQ. 8 and 9 can be used to update only the new objects. In addition, the 
sampling procedure ensures that the selected pairs contain at least one object from the 
incoming set. That is. two objects are selected at random so that at least one of these objects 
belongs to the incoming set. 

& Operation of the Present Invention 

25 The operation of the present invention with regard to visualizing and interactively 

processing chemical compounds in a non-linear map shall now be described with reference 
to a flowchart 302 shown in FIG. 3. Unless otherwise specified, interaction with users 
described below is achieved by operation of the user interface modules 108 (FIG. 1 ). 
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In step 304. the user selects one or more compounds to map in a new non-linear 
map. The user may select compounds to map by retrieving a list of compounds from a file, 
by manually typing in a list of compounds, and/or by using a graphical user interface (GUI) 
such as the structure browser shown in FIG. 5 (described below). The invention envisions 
5 other means for enabling the user to specify compounds to display in a non-linear map. For 
example, the user can also select compounds from an already existing compound 
visualization non-linear map (in one embodiment, the user drags and drops the compounds 
from the old compound visualization non-linear map to the new compound visualization non- 
linear map -- drag and drop operations according to the present invention are described 
10 below). 

In step 306, the user selects a method to be used for evaluating the molecular 
similarity or dissimilarity between the compounds selected in step 304. In an embodiment, 
the similarity/dissimilarity between the compounds selected in step 304 is determined (in 
step 308) based on a prescribed set of evaluation properties. As described above, evaluation 

15 properties can be any properties related to the structure, function, or identity of the 
compounds selected in step 304. Evaluation properties include, but are not limited to, 
structural properties, functional properties, chemical properties, physical properties, 
biological properties, etc., of the compounds selected in step 304. 

In an embodiment of the present invention, the selected evaluation properties may be 

20 scaled differently to reflect their relative importance in assessing the proximity (i.e., 
similarity or dissimilarity) between two compounds. Accordingly, also in step 306, the user 
selects a scale factor for each of the selected evaluation. Note that such selection of scale 
factors is optional. The user need not select a scale factor for each selected evaluation 
property. If the user does not select a scale factor for a given evaluation property, then that 

25 evaluation property is given a default scale factor, such as unity. 

Alternatively in step 306, the user can elect to retrieve similarity/dissimilarity values 
pertaining to the compounds selected in step 304 from a source, such as a database. These 
similarity/dissimilarity values in the database were previously generated. In another 
embodiment, the user in step 306 can elect to determine similarity/dissimilarity values using 

30 any well-known technique or procedure. 

In step 308, the map generating module 106 generates a new non-linear map. This 
new non- linear map includes a point for each of the compounds selected in step 304. Also, 
in this 
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new non-linear map, the distance between any two points is representative of their 
similarity/dissimilarity. The manner in which the map generating module 106 generates the 
new non-linear map shall now be further described with reference to a flowchart 402 in FIG. 4. 

In step 404, coordinates on the new non-linear map of points corresponding to the 
5 compounds selected in step 304 are initialized. 

In step 406, two of the compounds i, j selected in step 304 are selected for processing. 

In step 408, similarity/dissimilarity d i} between compounds i, j is determined based on 
the method selected by the user in step 306. 

In step 410, based on the similarity/dissimilarity d^ determined in step 408, coordinates 
1 0 of points corresponding to compounds i, j on the non-linear map are obtained. 

In step 412, training/learning parameters are updated. 

In step 414, a decision is made as to terminate or not terminate. If a decision is made 
to not terminate at this point, then control returns to step 406. Otherwise, step 416 is 
performed. 

15 In step 416, the non-linear map is output (i.e., generation of the non-linear map is 

complete). 

Details regarding the steps of flowchart 402 are discussed above. 
Referring again to FIG. 3, in step 312 the map viewer 1 12 displays the new non-linear 
map on an output device 1 16 (such as a computer graphics monitor). Examples of non-linear 
20 maps being displayed by the map viewer 1 12 are shown in FIGS. 6 and 7 (described below). 

In step 314, the user interface modules 108 enable operators to interactively analyze 
and process the compounds represented in the displayed non- linear map. These user interface 
functions of the present invention are described below. 

The present invention enables users to modify existing compound visualization non- 
25 linear maps (as used herein, the term "compound visualization non-linear map" refers to a 
rendered non-linear map). For example, users can add additional compounds to the map, 
remove compounds from the map, highlight compounds on the map, etc. In such cases, 
pertinent functional steps of flowchart 302 are repeated. For example, steps 304 (selecting 
compounds to map), 310 (generating the non-linear map), and 312 (displaying the map) are 
30 repeated when the user opts to add new compounds to an existing map. However, according 

to an embodiment of the invention, the map is incrementally refined and displayed in steps 310 
and 3 1 2 when adding 
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compounds to an existing compound visualization non-linear map (this incremental 
refinement is described above). 

9. User Interface of the Present Invention 

The user interface features of the present invention are described in this section. 
5 Various user interface modules and features are described below. Also, various 
functional/control threads (in the present context, a functional/control thread is a series of 
actions performed under the control of a user) employing these user interface modules and 
features are described below. It will be appreciated by persons skilled in the relevant art(s) 
that the user interface of the present invention is very flexible, varied, and diverse. An 

10 operator can employ the user interface of the present invention to perform a wide range of 
activities with respect to visualizing and interactively analyzing chemical compounds. 
Accordingly, it should be understood that the functional/control threads described herein are 
provided for illustrative purposes only. The invention is not limited to these 
functional/control threads. 

15 Preferably, the invention provides the following capabilities, features, and functions: 

displaying 2D and/or 3D chemical structures and/or chemical names; displaying compound 
collections and/or libraries; displaying components of structures (i.e. building blocks) of 
combinatorial libraries; visualization of compound collections and/or libraries as 2D and/or 
3D maps of colored objects. 

20 Also, the present invention allows the following: (1) browsing compound collections 

and/or libraries; (2) selection of individual compounds, collections of compounds and/or 
libraries of compounds; (3) selection of compounds generated in a combinatorial fashion via 
selection of their respective building blocks; (4) mapping, visualization, and/or linking of 
compounds onto and/or from 2D and/or 3D maps; (5) manipulation of the 2D and/or 3D 

25 maps such as rotation, resizing, translation, etc.; (6) manipulation of objects on the 2D and/or 
3D maps such as changing the appearance of objects (visibility, size, shape, color, etc.), 
changing position of objects on the map, and/or changing relationships between objects on 
the map; (7) interactive exploring of the 2D and/or 3D maps such as querying chemical 
structure, querying distance, selection of individual objects and/or areas of a map, etc. 
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Additional user interface features, functions, and capabilities of the present invention 
will be apparent to persons skilled in the relevant art(s) based on the discussion contained 
herein. 

As shown in FIG. 1, the invention includes a structure browser 1 10 and a map viewer 
5 1 12. At any given time, each of these can have multiple instances depending on the program 
use. 

9. 1 Structure Browser 

FIG. 5 illustrates a structure browser window 502 generated by the structure browser 
1 10. The structure browser window 502 includes a frame 504, a menu pane 506, and a group 
10 of labeled tabbed pages 508. Each tabbed page holds a molecular spreadsheet or a group of 
labeled tabbed pages. 

Each tab is associated with a compound collection (tabs 510) or a library, such as a 
combinatorial library (tabs 512). Selecting a collection tab 510 brings up a table of 
corresponding chemical structures. Selecting a library tab 512 brings up a group of tabbed 

15 pages corresponding to the sets of building blocks used to generate the library. Each of the 
library's tabbed pages works the same way as a compound collection tabbed page. In the 
example shown in FIG. 5, the tab 510 called "DDL0" is selected. DDL0 has three building 
block tabs 512, called "Cores," "Acids," and "Amines." The "Acids" collection tab is 
currently selected, so that a table 522 of the structures of the compounds in the "Acids" 

20 collection is shown. 

The browser window 502 includes a table 522, a slider 514, an input field 516, and 
two buttons: "Prev Page" 518 and "Next Page" 520. The slider 514, the input field 516, and 
the buttons 518, 520 facilitate browsing the content of the Acids table 522. If we consider the 
content of the table 522 as a contiguous ordered list of chemical structures (compounds or 

25 building blocks), that shown in the browser window 502 can be considered as a window 
positioned over the list. At any given moment this window displays part of the list depending 
on its position and the displayed part is equal to the size of the window, i.e., the number of 
cells in the table. Initially that window displays the top of the list. Moving the slider 514 
changes the position of the window over the list. Entering a value into the input field 516 

30 specifies the position of the window over the list. Pushing the "Next Page" button 520 moves 
the window one 
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window size down the list, pushing the u Prev Page" bunon 518 moves the window one window 
size up the list. 

The user can select compounds shown in the table 522 for various actions. For 
example, compounds can be selected using the browser window 502 as input for the generation 
5 of a new compound visualization non-linear map, or as input for adding compounds to an 
existing compound visualization non-linear map. Clicking with a left mouse button over a 
table cell selects or deselects the corresponding compound structure (toggling). Toggling on/off 
also changes the color of the cell, to indicate which cells have been selected. Selected 
structures are displayed on a first background color, and non-selected structures are displayed 
10 on a second background color. In the example of FIG. 5, certain cells 523 in table 522 have 
been selected. 

The menu pane 506 contains menus: File, Edit, Selection, Map, and/or other menus. 
The File menu facilitates file open/save, print, and exit operations. Edit menu contains 
commands for editing content of the table 522. The Selection menu provides options to 
15 select/deselect (clear) a current compound collection, a collection of building blocks of a 
combinatorial library, and/or all compounds. The Map menu includes commands for creating a 
map viewer and for displaying a selection of compounds in that map viewer. The latter option 
brings up a dialog window (FIG. 8), which allows the user to specify shape, color, and/or size 
of the selected objects, which will be used to represent the selected compounds on the map. 

20 9.2 Map Viewer 

A map viewer window 600 generated by the map viewer 1 12 is shown in FIG. 6. (also 

see FIGS. 6-10 and 13). A compound visualization non-linear map is displayed in a render 

area 614 of the map view window 600. 

In a preferred embodiment, the map viewer 112 is based on Open Inventor, a C++ 
25 library of objects and methods for interactive 3D graphics, publicly available from Silicon 

Graphics Inc. Open Inventor relies on OpenGL for fast and flexible rendering of 3D objects. 

Alternatively, the map viewer 112 can be based on a publicly available VRML viewer. 

Alternatively, any other software and/or hardware product allowing rendering of 3D 

objects/scenes can be used. 
30 In a preferred embodiment, 3D compound visualization maps of chemical compounds 

are implemented as Open Inventor 3D scene databases. Each map is build as an ordered 

collection of 
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nodes referred to as a scene graph. Each scene graph includes, but is not limited to, nodes 
representing cameras (points of view), light sources, 3D shapes, objects surface materials, 
and geometric transformations. Each chemical compound displayed on a map is associated 
with a 3D shape node, a material node and a geometric transformation node. 

5 Geometric transformation node reflects compound coordinates in the map. 3D shape 

node and material node determine shape, size and color of the visual object associated with 
the compound. Combinations of a particular shape, size and color are used to display 
compounds grouped by a certain criteria, thus allowing easy visual differentiation of different 
groups/sets of compounds. 3D shapes of the visual objects in the map include, but not limited 

10 to, point, cube, sphere, and cone. Color of a visual object in the map can be set to any 
combination of three basic colors: red, green and blue. Besides the color, material node can 
specify transparency and shininess of a visual object's surface. 

In an embodiment, an object's display properties (color, intensity of color, 
transparent, degree of transparency, shininess, degree of shininess, etc.) can represent 

15 physical, chemical, biological, and/or other properties of the corresponding compound, such 
as the cost of the compound, difficulty of synthesizing the compound, whether the compound 
is available in a compound repository, etc. For example, the larger the molecular weight of 
an object, the larger the size of the corresponding object in the display map. 

Each object or point displayed in the compound visualization non-linear map 

20 represents a chemical compound. Objects in the compound visualization non-linear map can 
be grouped into sets. 

By default, every time a set of compounds is mapped into a compound visualization 
non-linear map, a new set of graphical objects is created and added to the compound 
visualization non-linear map. All objects in a particular set can share the same attributes: 

25 shape, color, and size, thus providing an easy visual identification of the objects belonging to 
the same set or to different sets. 

A compound can be a member of several sets. In an embodiment, for a given 
compound, a different object is displayed in the compound visualization non-linear map for 
each set of which the compound is a member. In this case the objects in the compound 

30 visualization non-linear map that represent the compound as a member of each of the sets 
may overlap and only the 
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biggest object may be visible. In this case, a toggle sets feature (described below) may be 

used to reveal multiple set membership. 

The map viewer window 600 includes a frame 602, a menu pane 604, and a viewer 

module preferably implemented as an Open Inventor component (examiner viewer). The 
5 viewer module incorporates the following elements: (1) a render area 614 in which the 

compound visualization non-linear map is being displayed; (2) combinations of thumbwheels 

608, 610, 612, sliders, and/or viewer functions icons/buttons 620, 622, 624, 626, 628, 630, 

632; and (3) pop-up menus and dialogs 616, 702, 902 which provide access to all viewers 

functions, features and/or properties. 
10 The thumbwheels 608, 610 rotate the compound visualization non-linear map around 

a reference point of interest. Thumbwheel 610 rotates in the y direction, and thumbwheel 608 

rotates in the x direction. The origin of rotation (i.e., the camera position) is by default the 

geometric center of the compound visualization map 614 (render area), but can be placed 

anywhere in the compound visualization non-linear map. The compound visualization non- 
15 linear map can also be panned in the screen plane, as well as dollied in and out 

(forward/backward movement) via thumbwheel 612. 

The map view window 600 has several different modes or states, e.g. view, pick, 

panning, dolly, seek, and/or other. Each mode defines a different mouse cursor and how 

mouse events are interpreted. 
20 In the view mode, mouse motions are translated into rotations of the virtual trackball 

and corresponding rotations of the compound visualization non-linear map. The view mode 

is the default mode. 

In the panning mode, the compound visualization non-linear map is translated in the 
screen plane following the mouse movements. 
25 In the dolly mode, a scene is moved in and out of screen according to the vertical 

motions of the mouse. 

Seek mode allows the user to change the point of rotation (reference point) of a scene 
by attaching it to an object displayed in the compound visualization non-linear map. 

Pick mode is used for picking (querying) objects displayed in the compound 
30 visualization non-linear map. Picking an object in a 3D scene is achieved by projecting a 
conical ray from the camera through a point (defined by positioning and clicking the mouse) 
on the near plane of the 
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view volume. The first object in the scene intersecting with the ray cone is picked. As a 
response to a pick event (an object being picked by pressing the left mouse button over the 
object), a small window displaying the corresponding compound pops up while the left 
mouse button is pressed (see, for example, window 1302 in FIG. 13). The window will 
5 automatically disappear when the button is released. In order to keep the window on the 
screen, it is necessary to hold the shift key while releasing the mouse button. 

Switching between the above-described modes can be achieved by selecting a mode 
from a pop-up menu, by clicking on a shortcut icon/button, and/or by pressing and/or holding 
a combination of mouse buttons and/or keys on a keyboard. In a preferred embodiment, 
10 selecting a pointed arrow icon/button 620 switches to the pick mode. Selecting a hand 
icon/button 622 switches to the view mode; selecting a target icon/button 624 switches to the 
seek mode. Pressing and holding the middle mouse button switches to the panning mode. 
Pressing and holding the left and middle mouse buttons simultaneously switches to the dolly 
mode. 

1 5 Certain actions can be executed also via thumbwheels and/or sliders, e.g. turning the 

dolly thumbwheel 612 moves the scene in and out of the screen. Also, turning the X and/or 
Y rotation thumbwheels 608, 610 rotate the scene accordingly around the point of rotation. 

In a preferred embodiment, the right mouse button is reserved for the pop-up menus 
616, 902. Pressing the right mouse button anywhere over an empty* rendering area brings up 

20 the viewer pop-up menu 902. Pressing the right mouse button over an object brings up the 
object pop-up menu 616. 

The viewer pop-up menu 902 allows the user to select the mode (such modes are 
described above), change viewer properties (set up preferences, e.g. background color), 
toggle on/off sets of objects, and/or access any other viewer features. 

25 The object pop-up menu 616 allows the user to change an object's shape, color 

(material), and/or size, select the corresponding set of compounds, and/or define a 
neighborhood 3D area around the object (zoom feature, described below). In a preferred 
embodiment, ail changes made to an object automatically apply to ail other objects from the 
same set. The object's shape can be changed to one of the predefined basic shapes (e.g. dot, 

30 cube, sphere, cone). The object's material (color) is changed via a color dialog. The object's 
size is changed via a resize dialog. Any set of objects can be visible (toggled on) or hidden 
(toggled off). A toggle sets command 
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brings up a list of sets defined for the current map 640. Clicking on a set in the list 
(highlighting/clearing) toggles the set off and on. 

Invoking the zoom feature (via the pick neighbors command on the object pop-up 
menu 616, for example) creates a sphere 704 in the render area 614 (FIG. 7), which is 
5 centered on the object. The radius of the sphere 704 can be adjusted via a resize dialog 702 to 
select a desired neighborhood area around the object. AH objects (and corresponding 
compounds) encompassed by the sphere 704 are then selected, displayed in a different map, 
added to a new or existing set, dragged to a target (described below), and/or viewed in a 
structure browser window 502. 

10 The map viewer 112 is capable of maintaining an interactive selection of 

objects/compounds. All selected objects are visualized in the same shape, color, and/or size. 
In other words, selecting an object changes its shape, color, and/or size (e.g. to a purple 
cone), deselecting an object changes its shape, color and/or size back to the original 
attributes. Executing the select set command from the object pop-up menu 616 selects the 

15 whole set of objects this object belongs to. Alternatively, an individual object can be 
selected or deselected by clicking a middle mouse button over an object. The interactive 
selection of objects can be converted to a set of compounds and displayed in a structure 
browser window 502. The current selection can be converted into a set of compounds by 
invoking the save selection command from a selection menu, and/or it can be cleared by 

20 executing the clear selection command from the selection menu. 

9 J Interactivity of the Present Invention 

As should be apparent from the above, the present invention enables users to interact 
with the objects/compounds displayed in a compound visualization non-linear map. This 
interactivity provided by the present invention shall be further illustrated below. 

25 93 J Map Viewer as Target 

According to the present invention, a user can select a plurality of compounds from 
some source, and then add those compounds to a new or an existing compound visualization 
non-linear map being displayed in a map window 600. In this instance, the map window 600 
(or, 
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equivalently in this context, the map viewer 112) is acting as a target for an interactive user 
activity. 

This operation is conceptually shown in FIG. 14. A compound visualization non- 
linear map 1404 is being displayed in a map window 600. According to the present 
5 invention, the user can select compounds from a structure browser window 502, and then add 
those selected compounds (through, for example, well known drag and drop operations) to 
the compound visualization non-linear map 1404. Similarly, the user can select compounds 
from a compound database 122, or from a MS (mass spectrometry) viewer 1402, and then 
add those compounds to the compound visualization non-linear map 1404. 
10 According to an embodiment of the invention, new compounds are added to an 

existing compound visualization non-linear map by incremental refinement of the compound 
visualization non- linear map. Such incremental refinement is described above. 

9.3.2 Map Viewer as Source 

According to the present invention, a user can select a plurality of compounds from a 
1 5 map window 600, and then have those compounds processed by a target. In this instance, the 
map window 600 (or, equivalently in this context, the map viewer 1 12) is acting as a source 
for an interactive user activity. 

This operation is conceptually shown in FIG. 13. A user selects one or more 
compounds from the compound visualization non-linear map being displayed in the map 
20 window 600, and then drags and drops the selected compounds to a target. The described 
action is interpreted as a submission of the corresponding chemical structure(s) to the 
receiving target for processing. The receiving object can be anything that can handle a 
chemical structure: another map viewer 112, a structure viewer 110, a (molecular) 
spreadsheet 136, a database 120, an experiment planner 140, an active site docker 144, an 
25 NMR widget 130, an MS widget 134, a QSAR model 138, a property prediction program 
142, or any other suitable process. For example, dragging and dropping a compound onto an 
NMR widget would display this compound's NMR spectrum, either an experimental or a 
predicted one. 

The experiment planner is described in pending U.S. Patent Application titled 
30 "SYSTEM, METHOD AND COMPUTER PROGRAM PRODUCT FOR IDENTIFYING 
CHEMICAL 
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COMPOUNDS HAVING DESIRED PROPERTIES," Atty. Docket No. 1503.0200001, 
herein incorporated by reference in its entirety. 

The drag and drop concept described above provides a powerful enhancement of a 
3D mapping and visualization of compound collections and libraries. Any conceivable 
5 information about a set of chemical compounds can thus be easily accessed from the 
compound visualization non-linear map. For example, a map of compounds capable of 
binding to an active site of a given enzyme or receptor would benefit from the possibility to 
visualize how compounds from the different areas of the map bind to that enzyme or 
receptor. 

10 9.4 Multiple Maps 

According to the present invention, it is possible to create multiple visual maps for 
any given set of collections and/or libraries of chemical compounds. Multiple visual maps 
can be based on the same and/or different non-linear maps. Visual maps based on the same 
non-linear map can display different subsets of compounds and/or present different views of 

15 the same set of compounds (e.g. one visual map can display an XY plane view and another 
visual map can display an orthogonal, YZ plane view). Visual maps based on different non- 
linear maps can visualize the same set of compounds on different projections, for example, 
maps derived from different similarity relations between these compounds. 

If a compound is mapped on multiple visual maps, the visual objects representing the 

20 compound on the different maps can be crossl inked. Crosslinking means that any 
modifications made to a visual object in one of the visual maps will be automatically 
reflected into the other visual maps. For example, if an object is selected on one of the visual 
maps, it will be displayed as selected on the other visual maps as well. In fact, all objects on 
all maps can be crosslinked provided that they represent the same chemical compounds. 

25 Multiple visual maps can be also crosslinked in a way that mapping any additional 
compounds onto one of the visual maps will automatically map the same compounds onto the 
crosslinked maps. 

10. Examples 
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The present invention is useful for visualizing and interactively processing any 
chemical entities including but not limited to small molecules, polymers, peptides, proteins, 
etc. It may also be used to display different similarity relationships between these 
compounds. 

The present invention has been described above with the aid of functional building 
blocks illustrating the performance of specified functions and relationships thereof. The 
boundaries of these functional building blocks have been arbitrarily defined herein for the 
convenience of the description. Alternate boundaries can be defined so long as the specified 
functions and relationships thereof are appropriately performed. Any such alternate 
boundaries are thus within the scope and spirit of the claimed invention. These functional 
building blocks may be implemented by discrete components, application specific integrated 
circuits, processors executing appropriate software and the like or any combination thereof. 
It is well within the scope of one skilled in the relevant art(s) to develop the appropriate 
circuitry and /or software to implement these functional building blocks. 

While various embodiments of the present invention have been described above, it 
should be understood that they have been presented by way of example only, and not 
limitation. Thus, the breadth and scope of the present invention should not be limited by any 
of the above-described exemplary embodiments, but should be defined only in accordance 
with the following claims and their equivalents. 
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What Is Claimed Is: 

1. A computer implemented method for visualizing data relating to chemical 
compounds, comprising the steps of: 

( 1 ) enabling a user to select a plurality of compounds to map; 

(2) enabling the user to select a method for evaluating the 
similarity/dissimilarity between said selected compounds; 

(3) generating a non- linear map in accordance with said selected compounds and 
said selected method, said non-linear map having a point for each of said selected 
compounds, wherein a distance between any two points is representative of 
similarity/dissimilarity between corresponding compounds; and 

(4) displaying at least a portion of said non-linear map. 

2. The method of claim 1, further comprising the step of: 

(5) enabling the user to interactively analyze compounds represented in said 
non-linear map. 
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